Velocity-modulation control of electron-wave propagation in graphene 
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Wave propagation control by spatial modulation of velocity has a long history in optics and 
acoustics. We address velocity-modulation control of electron wave propagation in graphene and 
other two-dimensional Dirac-electron systems, pointing out a key distinction of the Dirac-wave case. 
We also propose a strategy for pattern transfer from a remote metallic layer that is based on many- 
body velocity renormalization. 
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Introduction — Control of electromagnetic and mechani- 
cal wave propagation by spatial modulation of wave ve- 
locity [U [2] has been studied for many decades in op- 
tics and acoustics, originally in relatively simple multi- 
layer structures which have a wealth of practical ap- 
plications and more recently in sophisticated two- and 
three-dimensional photonic and phononic crystal struc- 
tures which can have gaps between transmission bands. 

Recent advances [3] in the isolation and control of sin- 
gle and few-layer graphene electron systems motivate a 
close examination of velocity-modulation control in this 
material. Graphene is an allot rope of carbon atoms 
tightly packed in a two-dimensional (2D) honeycomb lat- 
tice. At energies near the Fermi energy of a neutral sys- 
tem, electron waves in graphene are described by a 2D 
massless Dirac equation and, like electromagnetic and 
mechanical waves, travel with a velocity that is indepen- 
dent of wavelength. Because of this property, the analo- 
gies of electron transport with both optics [4] and acous- 
tics are stronger than in the conventional non-relativistic 
electron- wave case [5]. In this article we first consider 
the propagation of massless Dirac fermion (MDF) waves 
through a medium with a position-dependent velocity, 
highlighting a key distinction between a Dirac wave and 
electromagnetic or acoustic waves. We then discuss a 
non-invasive strategy for achieving velocity modulation 
in graphene without any direct physical contact to the 
sample, by transferring a spatial pattern from remote 
metal layers via many-body velocity renormalization. 
Scattering of MDFs against a velocity harrier — The in- 
fluence of velocity variation on propagation is best illus- 
trated by the simplest example, transmission through a 
velocity harrier [6] , as illustrated in the inset in Fig. [l] 
We flrst solve this scattering problem, highlighting a key 
distinction between a Dirac wave and electromagnetic or 
acoustic waves. We consider MDFs in a medium in which 
the Fermi velocity v changes as a function of the 2D po- 
sition r: V = v{r). The massless Dirac equation in this 
case reads [7 

n^f{r) = -ih^v{r) a Vr [^/v{r)^f {r)] = E ^(r) , (1) 



where ^(r) = (^a(^), ^b(^))^ is a two-component 
spinor, ^a(^) and ^b(^) are the honeycomb sublattice 
components of the electron wave, and cr = {cFx^cFy) is a 
2D Pauli matrix vector. In using the Dirac equation we 
are assuming that velocity variations are slow on a lat- 
tice constant scale. In this limit spin and valley degrees 
of freedom play a passive role. Note that the Hamiltonian 
in Eq. ([l]) does not differ from its uniform system coun- 
terpart TYmdf = —ihva • Vr merely by the replacement 
V v{r): as pointed out by Peres in Ref. |7| this pre- 
scription would lead to a non-Hermit ian operator. It is 
nevertheless convenient to introduce the auxiliary spinor 
^{r) = y/v{r)'^{r) which satisfles 

- ihv{r)a ■ Vr^{r) = E ^(r) . (2) 

For the barrier problem illustrated in Fig. = v{x) 

changes only along the x direction and momentum along 
the y direction is a good quantum number. 
We solve Eq. ([2| for a simple velocity barrier 

{^1, if X < 
^2, if < X < D . (3) 
vi^ \i X > D 

This step-like model is justifled when the Fermi wave- 
length A is much larger than the characteristic width 
over which v{x) changes, as discussed at length for 
the potential-barrier case in Ref. |8l For a given y- 
direction wavevector ky we are left with two coupled one- 
dimensional flrst-order differential equations: 

- ihv{x) T 0A(B)(^) = E (t)B{A){x) . (4) 

where ^{r) = (j)[x)e^^yy and the spinor = 
((/>a(^), (/>b(^))^- The flrst order equations for the spinor 
components can [because v{x) is piecewise constant] be 
combined into a second order equation satisfled by both: 





" E ' 






Jiv{x)_ 





2 



1.0 



0.8 



0.6 



0.4 



0.2 



0.0 









\\ :\ 


' / \ ' 

/ \ 








/ \ 


y 

Vi 
















I 

I 

V \ 

\ 1 






X j 








= 0.2 i 








= 0.5 i 






e 


= 1.2 i 






— e 


= 1.5 ' 














\, 1 ^Oii 



7r/4 

ei 



37r/8 



7r/2 



FIG. 1: (Color online) Inset: Cartoon of a velocity barrier. 
The velocity of the massless carriers changes along the x di- 
rection according to the simple functional form defined by 
Eq. (|3|. Main panel: Transmission probability T^,^(^i) as a 
function of angle of incidence Oi for velocity ratio ^ = V2/V1 
equal to 0.2, 0.5, 1.2, and 1.5. The angle of incidence is defined 
so that ^1=0 corresponds to propagation perpendicular to 
the interface. These data have been obtained for kD — rj = An 
where k is the incoming wavevector and D is the barrier thick- 
ness. This value is typical 8 and corresponds for example to 
k — 27r/A with A = 50 nm and D — 100 nm. Note that for 
^ > 1, the transmission probability drops rapidly toward zero 
when the critical incidence angles, indicated by vertical bars, 
are exceeded, (^c ~ 0.99 rad for ^ = 1.2 and ~ 0.73 rad for 
e = 1.5.) 



The solutions of Eq. ([5| can be written as 

(^^ik^x rj^^-ik^x^^^ if X < 



(t)A{x) = <^ (ae*^-^ + 6e-*^-^), if < x < D , (6) 
[ te^^-^, \ix> D 

where the constants a, 6, r, and t are to be determined, 



knr; 



Qx 




(7) 



sin^((9i) 



k = E / ihvi)^ = V2/vi\^ the velocity ratio, and Oi is the 
angle of incidence, i.e. tan(^i) = ky/kx- The correspond- 
ing expression for ^b(^) can be obtained from Eq. (|6| us- 
ing Eq. Q. InEq. Q we can identify O2 = tan ^{ky/qx) 
as the angle of refraction and thereby obtain a quantum 
version of the famous Snell- Descartes law of geometrical 
optics: 



sin(6>i) 
sin(^2) 



V2 



(8) 



For ^ > 1 and Oi > 6^ = sin ^(1/^), the classical total 
internal reflection angle, Qx is imaginary and we expect 



negligible transmission through thick barriers. In this 
case the classical correspondence fails and the refraction 
angle is not well defined. 

Explicit evaluation of the four coefficients in Eq. (|6| 
requires matching conditions at the two interfaces which 
we obtain by the following argument. Dividing both sides 
of Eq. Q by v{x) and integrating across either inter- 
face implies that the auxiliary spinor <l> is continuous, 
and therefore that the physical ^ satisfies the following 
matching conditions: 



*(0+,y) = -^ *(0-,y) 



(9) 



These discontinuities in ^ guarantee that the divergence 
of the local current J(r) = v{r)^\r)a^ {r) vanishes. 
Using these matching conditions we are able to obtain 
an explicit expression for the transmission probability 



cos^ {ei)[l-e sin^ (^1)] 

Cr,,5(6'l) 



(10) 



cos2 (6'i)[l-e2 sin^ (6*1)] 



where rj = kD and Cr, 5(^1) 
(l-0'sin2 (g,7?A)sin2 {9,). 

In Fig. [1] we plot the transmission probability T^^^ as 
a function of Oi for = 47r at several different values 
of the velocity ratio ^. Note that a velocity barrier is 
always perfectly transparent, T^^^ = 1, for normal inci- 
dence (^1 = 0) as in the standard Klein problem [8]. This 
property of Dirac-wave propagation through a velocity 
barrier establishes a qualitative difference between the 
present case and the familiar electromagnetic and acous- 
tic cases, and opens up new ground for the invention of 
spatial patterns with desirable transmission properties. 
It is also an important addition to the obvious differ- 
ence in velocity- wavelength relationship in distinguishing 
Dirac-wave propagation from Schrodinger-wave propaga- 
tion. In Fig. [2] we plot the integrated transmission, 

T,(o=- r^'d^iT,,^(^i), (11) 

as a function of the velocity ratio ^. A sharp change in 
behavior is visible at ^ = 1 which foreshadows the total 
internal reflection properties of classical waves. Indeed, 
in the 77 ^ 00 limit it is easy to prove that for ^ = 1+ 
T^(0 ^ 2i9c/7r ^ 1 - 2y2vT^/7r, which is a non- 
analytic function of ^. For this reason even a slight mis- 
match in velocities can produce a large electron transport 
signal. 

Molding electron flow — Doped or gated graphene sheets 
are normal Fermi liquids [TOl |Tl] but have a number 
of unusual quantitative features p!2^ in their correlation 
physics which might provide an attractive route to ve- 
locity modulation, as we now explain. We first con- 
sider a grounded metal plane placed close to a graphene 
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FIG. 2: (Color online) The integrated transmission as 
a function of the velocity ratio ^ for rf — tt [solid (black) 
line] and 77 = 47r [dashed (red) line]. The dotted (blue) 
line shows the classical limit Too{^ 1^). The dash- 
dotted (green) line shows the integrated transmission for 
a linearly-polarized (along the ^-direction) electromagnetic 
wave scattering against a non-magnetic barrier (in this case 
^ = V^ei/e2, where a are the dielectric constants of the media 
composing the barrier, and 77 = ijjD /c). 

sheet. The presence of the metal does not directly shift 
the chemical potential, like a biased gate would, but be- 
cause electron-electron interactions between MDFs are 
expected to be substantially screened and thus weaker, 
has a rather large impact on the renormalized Fermi 
velocity of the quasiparticles close to the Fermi en- 
ergy that are important in transport. We show below 
that quasiparticles under the screening plane move at a 
speed v'^ that is smaller than in an isolated graphene 
sheet [Ml [15]. Using a single ground metal plane located 
close to a graphene sheet will thus lead to a very simple 
realization of the velocity barrier illustrated in the inset 
in Fig. [1] with ^ < 1. To realize a velocity barrier with 
^ > 1 one instead needs to use two metal gates located 
on top of the regions with x < and x > D. The area on 
top of the strip < x < D must instead be left empty. 
In order for our abrupt-interface velocity-barrier calcu- 
lation to be relevant, the distance to the metal gate d 
would have to be smaller than the Fermi wavelength A, 
but transmission properties will be similar even if this 
condition is not satisfied. Any shape of velocity modu- 
lation can be achieved by transferring a suitable spatial 
pattern from the remote lithographically-designed metal 
layer. 

The effect of a metal gate on the quasiparticle veloc- 
ity can be estimated quantitatively by evaluating the 
quasiparticle self-energy S of an interacting MDF sys- 
tem near the quasiparticle pole. We have generalized 
the GW theory [16 calculations described in Ref. tl4| 



replacing the bare Coulomb interaction with the corre- 
sponding expression appropriate for a 2D electron sys- 
tem close to a perfect metallic screening plane: Vd{q) = 
27re2[l -exp(-2(7d)]/g [17 . 

In Fig. ^) we report numerical results for v'^/v as a 
function of electron density n for several values of and 
in Fig. [s})) as a function of d for fixed density. We see 
that a substantial velocity contrast can be induced by 
metallic gates that are tens of nm's from the graphene 
plane. The effect of the gate can extend much further 
if it is separated from the graphene by a dielectric with 
er > 1. 

In addition to the velocity modulation the remote 
metallic layer will produce a shift in chemical potential. 
The reason lies in the fact that the Fermi energy quasi- 
particles whose transmission properties we study here 
satisfy, strictly speaking, a Dyson equation [14^, not a 
single-particle massless Dirac equation. The chemical 
potential shift acts on the quasiparticles exactly like an 
external potential acts on a free particle. These shifts are 
small however, and can be compensated by biasing the 
patterned metallic plane although it might be difficult to 
completely separate velocity patterning from potential 
patterning in real experiments. 

Discussion — Most of the considerations outlined in this 
paper apply equally well to any system in which mat- 
ter waves satisfy a massless Dirac equation, for example 
to the surface states of topological insulators Re- 
cently it has been proposed [19 i [20l [21] that MDFs can 
also be realized in any standard 2D electron gas (2DEG), 
when appropriately nanopatterned. Similar proposals to 
realize the Dirac spectrum have been discussed in the 
contexts of ultracold atoms in optical lattices [22] and 
photonic crystals [23]. Velocity modulation can be re- 
alized in these systems as well: a miniband structure 
is imprinted on a 2DEG subjected to a long- wavelength 
periodic external potential {i.e. a lateral superlattice) 
with hexagonal symmetry. If suitable conditions are sat- 
isfied [20] , isolated Dirac points described by simple MDF 
Hamiltonians can appear in this miniband structure. As 
shown in Ref. 20 the Fermi velocity in these systems is 
quite sensitive to the strength of the periodic potential. 
By patterning the surface of a 2DEG in such a way to 
create three regions along a given direction in which the 
strength of the external periodic potential changes one 
can achieve a velocity barrier similar to the one sketched 
in the inset in Fig. [l] 

In summary, we have calculated the transmission prob- 
ability of massless Dirac fermions through a model "ve- 
locity barrier" and showed how electrons flowing through 
it obey the Snell-Descartes law of optics. We have also 
discussed a practical strategy for achieving substantial 
velocity modulation without damaging the graphene by 
exploiting the influence of a remote metallic layer on 
many-body renormalization of the quasiparticle velocity. 
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FIG. 3: (Color online) Panel a) Renormalized quasiparticle 
velocity (in units of the bare velocity v) in graphene as 
a function of the electron density n (in units of 10^^ cm~^) 
in the presence of a ground perfect-metal plane located at a 
distance d from the sheet (see inset). The (red) filled circles 
refer to an isolated sheet (in the absence of the perfect-metal 
plane). The other data refer to finite and increasingly larger 
(from top to bottom) values of d. Note how jv is largely 
affected by the presence of the ground metal plane. Panel b) 
shows v'*' jv as a function of d (in nm) for different values of 
n\ V* /v is largely suppressed even when the gate is quite far 
away from the graphene sheet. 
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